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Abstract. We present the numerical simulations of co- 
herent inverse Compton scattering (CICS) in a highly 
magnetized plasma process by means of a full three di- 
mensional particle in cell code (PIC), which is mass and 
energy conservative. We used the parameters of a pulsar 
magnetosphere where CICS is one of the most promising 
models for the generation of the observed highly coherent 
radio emission. First we show details of the onset and time 
evolution of strong Langmuir turbulence driven by a rela- 
tivistic electron beam penetrating a strongly magnetized 
background plasma. The Langmuir turbulence acts as self- 
generated wiggler fields which bunch the beam electrons 
thereby inducing strong coherent emission of the bunches 
at frequency 7^ times the plasma frequency. The emitted 
power is about 10 GW in a few nanoseconds. This radia- 
tion is interpreted in terms of inverse Compton scattering 
on nonlinear density fluctuations. CICS is the longitudi- 
nal version of a laboratory free electron laser and is appli- 
cable in strongly magnetized plasmas like pulsars. 



1. Introduction 



to nonlinear amplitudes (Benford and Weatherall 1991) 



(Kato et al. 1983), (WeatheraU 1988). The interaction of 



the REB with the electrostatic fluctuations of the Lang- 
muir waves results in very strong coherent electromagnetic 

In the rest 



(Kato et al. 1983), (Levron ct al. 1987) 



frame of the beam electrons, they experience the Doppler 
shifted electrostatic wave moving towards them. The non- 
linear waves force the REB to wiggle and to emit dipole 
radiation of the same frequency as the waves. In the labo- 
ratory frame this Hertz' dipole radiation is again Doppler 
boosted and strongly beamed in forward direction due to 
the relativistic lighthouse effect. This interaction can be 
viewed as coherent inverse Compton scattering (C ICS) 
off Langmuir waves (Benford and Weatherall 1991) , Co- 



herence occurs due to the bunch structure of the beam 
which means a phase coupled electrostatic field and den- 
sity modulation. 



2. Computational Details 



Since CICS 
candidate for 
ent radio emission from pulsars (Kunzl et al. 1998) 



has been discussed 
the origin of the 



A strong source for coherent radiation in laboratory is (Melrose and Gedalin 1999) , (Melrose 2000) , we use the 



the free electron laser (FEL) (O'Shea and Freund 2001) 



to be promising 
extremely coher- 



Its principle is the interaction of a relativistic electron 
beam (REB) and a magnetic wiggler system. The beam 
of high energy electrons with energy 7TOoC^ passes be- 
tween permanent magnets of alternating polarity and pe- 
riodicity Aw- The energy of the REB is efficiently con- 
verted by electron bunching. In the observer's frame the 
bunches emit the relativistically Doppler-shifted wave- 
length A = Aw /t^ ■ A comparable effect can be driven by a 
REB penetrating a strongly magnetized plasma. In phase 
space such a configuration corresponds to a two-stream- 
instability. In the plasma version the free energy of the 
emitted radiation is supplied by the REB, as it is in the 
case of the FEL, but the positive gradient of the momen- 
tum distribution function dj /dp > is the source of the 
instability which excites strong density fluctuations. In a 
strongly magnetized plasma such density oscillations are 
predominately longitudinal Langmuir waves which grow 



physical parameters of a pulsar magnetosphere as an in- 
structive application. Radio observations indicate that at 
about 500 km from the strongly magnetized neutron star 
the coherent pulsed radio emission is produced. The co- 
herence is proved by pulse substructures on time scales 
down to a few tens of nanoseconds with increasing flux 
densities (Rickett 1995), In these regions of a pulsar mag- 
netosphere typical values are: background electron den- 
sity rie = 10^^ m"'^, beam electron density equal to the 
plasma density nb = 10^^ m"'^, beam electron Lorenz fac- 
tor 7 = (7/3 = 2) and an extremely strong background 
magnetic field B of about 1000 T. Especially the strong 
guiding field is important for the efficiency of the CICS- 
process. It reduces the electron dynamics to the spatial 
dimension along the magnetic field lines which greatly im- 
proves the stability of the excited Langmuir waves and 
thus leads to strong a nd stable electrostatic wiggler fields 
(Pelletier et al. 1988)| . For the temperature of both the 
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Fig. 1. The initial configuration of electron phase space 
including the beam. The initial disturbance of the back- 
ground plasma is heavily exaggerated for reasons of better 
display. 



duce any finite charge or current densities, thus leaving 
the background plasma initially homogeneous, neutral and 
current free. 

The boundary conditions during the simulation are 
quite simple. In the transversal directions x and y we chose 
periodic boundaries and in the longitudinal direction z the 
box is open, which means that plasma particles and fields 
can leave the box unaffected. In addition at the lower z 
plane [z = 0) an electron component is continuously in- 
jected with the density profile mentioned above, a mean 
7/9 of 2 and the temperature Tb. 

The simulation runs for 600 time steps with 5t ~ 
2.41ns. The total time simulated was 1.44 /xs, after which 
the growing nonlinear Langmuir waves saturate and a 
quasi stationary situation had developed as it is shown 
in the next section. 



beam- and background electrons values we take the value 
of Tc = Tb = 100 eV, typical polar cap temperature of 
neutron stars derived from their thermal X-ray-emission 



(Pavlov and Petekhin 1995) The PIC-code used for this 



numerical simulation is fully three dimensional and it con- 
serves mass and energy. In the simulation the direction 
of electron beam propagation 7/3 and of the magnetic 
field B is the positive z~direction, which is also called 
the longitudinal direction. The simulated box has a nu- 
merical extension of 80 x 80 x 200 grid points and a phys- 
ical extension of 100 m x 100 m x 250 m, which corre- 
sponds to the size of the expected features. The resolu- 
tion of Sx = Sy = Sz = 1.25 m is more than sufficient 
to resolve the expected maximum wavelengths of approx- 
imately Aftad = 27r/fcRad = 27rc/7^ijjpc ~ 6.7m. A particle 
density of 10^^ m"'^ is approximated by 10 quasi particles 
per cell, leading to 16 million background quasi electrons 
that are accompanied by a varying number of beam quasi 
electrons of not more than 2.5 million, which eventually 
enter the computational box. Those numbers are doubled, 
since for quasi neutrality an equally large number of pos- 
itive particles are required. 

The initial condition is given by a homogeneous back- 
ground plasma of density rie, temperature Tg and a longi- 
tudinal magnetic field i?z = 0.1 T. This value is found to 
be sufficiently high to represent a much stronger field of 
1000 T. Initially there is no beam in the computational 
box. The beam is injected later at the z = plane of the 
box with a density profile of n (r) = nb/ (l -I- exp [^^^] ) ■ 
Here r denotes the distance from the x~y center of the 
box and i? = 20 m and AR = 3 m give the radius and the 
sharpness of the electron beam. The homogeneous back- 
ground plasma is disturbed in phase space (see Fig. |l]) 
in order to accelerate the excitation of Langmuir waves 
by giving the system a hint of the relevant length scales. 
We used the wave length of the fastest growing Langmuir 
mode, which follows from an analytic calculation of the 
linear growth rate. The disturbance is chosen not to in- 



3. Results 

The beam particles entering the computational box inter- 
act with the faint traces of fluctuations we introduced in 
our initial setup and excite the appropriate modes to grow. 
In Fig. ||on page |^ the temporal evolution of the electrons 
phase space in z-direction is presented by showing phase 
diagrams for the time steps 150, 300, 325 and 400. In the 
first plot (upper left) the penetrating electron beam can 
be seen together with the background plasma, both are 
still nearly undisturbed. In the following plots the develop- 
ment of strong Langmuir waves can be observed. After 781 
ns, the disturbance in the beam and the background has 
grown so strong, that both components start to mix, i.e. 
there is no distinction between beam and background elec- 
trons anymore. Beam electrons are strongly decelerated 
even until u ~ 0, whereas background electrons are accel- 
erated up to a Lorenz factor of 5. This can be seen best 
in Fig. on page'|[ where each of the both components 
are plotted in separate diagrams for the time step 400. In 
Fig. left the phase space of the beam electrons is shown. 
Clearly visible are the ring structures of fully developed 
Langmuir waves. Some of beam electrons are accelerated 
up to a 7/9 of 4-5 as it is the case for the background elec- 
trons, which is shown in Fig. |^ right. On the other hand 
the beam electrons are also decelerated to zero momen- 
tum during their interaction with the Langmuir waves. 
The extremely strong accelerations are due to the strong 
electrostatic fields of the excited Langmuir waves. Since 
the Langmuir turbulence structure in Fig. ^ moves with 
a much lower speed (which is even continuously decreas- 
ing throughout the simulation) than the yet unaffected 
beam electrons, it acts as the wiggler field that stimu- 
lates the CI CS-emission. The electrostatic fiuctuations 
directly correspond to electron density fluctuations (solid 
curve in Fig. |^). which reach a value of ~ -1-100%/ — 50% 
oriented at n ~ ric + — 2 ■ 10^^ m""^. The structure of 
the electron density is better seen in Fig. || (page ^ and 
Fig. H (page 0), where the temporal evolution of contour 
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Fig. 2. The electrostatic mode energy density is plotted 
against wavenumber k and time t. 



surfaces at n = 1.5 • 10^^ m"^ and rie = 1.9 
of a cut through the electron beam is shown. The electron 
density varies between 4 and 1 times 10^^ m""^, which im- 
plies the strong electrostatic fields. It can be clearly seen, 
that the beam decays in a series of " pancakes" . Such pan- 
cakes are absolutely necessary to explain the coherent na- 
ture of the pulsed radio emission of neutron stars. Why? 
Since the emission from relativistic particles is confined to 
a forward cone with half angle ^ 1/j and the emission 
is nearly along the magnetic field. Thus, in Fourier space 
bunching emission corresponds to the component k± to 
the fields being smaller than the component A:|| along the 
field line by a factor I/7. This corresponds to flat pan- 
cake shaped bunch with the normal within an angle I/7 



of B (Melrose 2000), This is exactly what we observe in 



our simulations! We briefly note that our simulations are 
the first ones which show this kind of pancake bunches! 
The number of particles in such a pancake corresponds to 
the the maximum coherence achievable for CICS and it 
should be noted, that in our simulation the thickness of 
such a pancake is significantly larger (> 1 m) than the De- 
bye length which is usually used in analytical models. A 
typical volume of one bunch is about lOm^ giving a total 
number of particles per bunch of about lO^""*. Our simu- 
lations prove, that a much higher coherence is achievable 
than expected from analytical models. 

The growth of Langmuir waves can be deduced by the 
Fourier transform of the electrostatic field E^- In Fig. || 
the energy density of the electrostatic modes is plotted 
versus the corresponding wavenumber k and time t. The 
k ~ 0.25 m"^ mode is growing by four orders of magnitudes 
within some ten nanoseconds. It should be noted however, 
that the fastest growing mode is slightly off compared to 
the initially given wave number can be seen in Fig. ^ as the 
small peak at early times. This means, that for our sim- 
ulation the analytical value for the fastest growing mode 
is not precisely correct and it proves, that although the 
initial disturbance has given a hint to the system it has 




Fig. 3. Structure of the emitted pointing flux as 3D con- 
tour surface of the value 2.5 • 10* W/ m^ (left) and as a cut 
through the center of the beam at a; = 50m (right). 
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Fig. 4. Dispersion relation (left) and emission character- 
istic (right) of the emitted radiation. For the left plot the 
angle 6m, the angle of maximum intensity, has been used. 
The curve in the left plot gives the dispersion branch of 
electromagnetic radiation in a plasma. The right diagram 
presents the intensity of the radiation for a given in 
a polar plot. 9 is the angle between beam direction and 
direction of emission. The values in both diagrams are 
normalized to a value of one. 



not forced it towards an unphysical solution. We note that 
the wave energy density is distributed to higher k when 
it exceeds a value of about 10~^ Jm^, driven by nonlinear 
wave-wave interactions. 

It is important to recognize that the generated elec- 
trostatic fluctuations react on the beam electrons sim- 
ilar to the wiggler fields in the laboratory FEL, with 
one major difference. In our case the wave wiggling self- 
consistently excited is along the propagation direction of 
the beam, whereas in the FEL-case the external mag- 
netic wiggling is perpendicular to the beam propaga- 
tion. The strongest emission originates in regions with 
the largest density gradients Fig. |[ Every de- and ac- 
celeration region forces the beam electrons to emit Hertz' 
dipole radiation, which is beamed in forward direction, 
as can be seen in Fig. ^ (right). 9m represents the angle 
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Fig. 5. Total integrated poynting flux leaving the com- 
putational box at the top (z ~ 250 m) against time. The 
+ symbols show the total power injected as kinetic energy 
of the beam electrons. 

of maximum intensity, expected from theoretical consid- 
eration (tan^M = 1/27/9) with a 7/3 = 2. The striking 
resemblance with a relativistically beamed Hertz' dipole 
is obvious, but the emission characteristic in our simula- 
tion also gives precisely the expected opening angle of the 
emission cone, which is proposed for a CICS-process. In 
Fig. ^ (left) the "spectrum" of the emitted radiation is 
shown, which has a strong peak at w = Swpe and a couple 
of minor peaks at 4wpe and Swpe due to nonlinear wave 
coupling. Again we find rather precisely the theoretically 
expected value of w = 7^a;po. The total power emitted 
by the simulated computational box is shown in Fig. 
It is the escaping poynting flux at the top of the compu- 
tational box (z = 250 m) integrated over the top surface. 
After about 850 ns the emitted power rises within some 
ten nanoseconds by more than 6 orders of magnitude to 
a value of 1 GW. Afterwards the total output power in- 
creases more slowly by one more order of magnitude and 
reaches a maximum of 10 GW continuous radiative power 
at the end of the simulation. This is a significant fraction 
of the power injected by kinetic energy of the beam par- 
ticles, shown in Fig. ^ as -|- symbols at a constant value 
of 50 GW. The rest of the energy leaves the box as kinetic 
energy of beam electrons, which would probably be also 
emitted if the computational would have a greater size in 
z-direction. 



electric field the inverse Compton scattering is coherent. 
Such a process is a very promising candidate for the co- 
herent radio emission of pulsars. 
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4. Discussion 

We have shown that in a strongly magnetized plasma a rel- 
ativistic electron beam can be forced to emit highly coher- 
ent radio emission by self-induced nonlinear density fluc- 
tuations. Such slowly moving nonlinear structures oscillate 
with the local plasma frequency at which the relativistic 
electrons are scattered. Beam electrons dissipate a signif- 
icant amount of their kinetic energy by inverse Compton 
radiation at a frequency of about ^'^ujpc. Since the beam is 
sliced into pancake structures which experience the same 
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Fig. 6. Electron phase space in ^-direction after 150, 300, 325 and 400 timesteps. The quasi particles are plotted as 
dots according to their corresponding position in z (ordinate) and their momentum pz (left abscissa) . The solid curves 
present the particle density of the electrons (right abscissa). 
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Fig. 7. Electron phase space in z-direction after 400 timesteps separated for beam (left) and background (right) 
electrons. 
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Fig. 8. Three dimensional contour plots of the electron particle density. For the timesteps 300, 400 and 500 (left to 
right) the surfaces are shown, where the density is n = 1.5 • 10^^ m"'^ (top) and Uc = 1.9 • 10^^ m"'^ (bottom) respectively. 
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Fig. 9. two dimensional contour plots of the electron particle density at timesteps 300,350,400 and 500. Shown is a 
cut through the center of the penetrating electron beam at x = 50 m. 



